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ABSTRACT 

We use a highly homogeneous set of data from 132 early-type galaxies in the Virgo and Fornax clus- 
ters in order to study the properties of the globular cluster luminosity function (GCLF). The globular 
cluster system of each galaxy was studied using a maximum likelihood approach to model the intrin- 
sic GCLF after accounting for contamination and completeness effects. The results presented here 
update our Virgo measurements and confirm our previous results showing a tight correlation between 
the dispersion of the GCLF and the absolute magnitude of the parent galaxy. Regarding the use of 
the GCLF as a standard candle, we have found that the relative distance modulus between the Virgo 
and Fornax clusters is systematically lower than the one derived by other distance estimators, and 
in particular it is 0.22mag lower than the value derived from surface brightness fluctuation measure- 
ments performed on the same data. From numerical simulations aimed at reproducing the observed 
dispersion of the value of the turnover magnitude in each galaxy cluster we estimate an intrinsic dis- 
persion on this parameter of 0.21mag and O.lSmag for Virgo and Fornax respectively. All in all, our 
study shows that the GCLF properties vary systematically with galaxy mass showing no evidence for 
a dichotomy between giant and dwarf early-type galaxies. These properties may be influenced by the 
cluster environment as suggested by cosmological simulations. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: star clusters — 
globular clusters: general 



1. INTRODUCTION 

The distribution of globular cluster (GC) magnitudes 
has the remarkable property that it is observed to peak 
at a value of My ~ —7.5 mag in a near universal fashion 
(e.g., Jacoby et al. 1992, Harris 2001, Brodie & Strader 
2006). This distribution, usually referred to as the GC 
luminosity function (GCLF), has been historically de- 
scribed by a Gaussian. By virtue of its near universality, 
the derived mean or "turnover" magnitude fi has seen 
widespread use as a distance indicator (e.g. Seeker 1992, 
Sandage & Tammann 1995), even though some disper- 
sion and discrepant results have been reported in the 
literature (see discussion in Ferrarese et al. 2000a). 

There is nevertheless no solid theoretical explanation 
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for the observed universality of the turnover magnitude. 
The luminosity function is a reflection of the more fun- 
damental mass spectrum of the GCs, and as such the 
"universal" turnover magnitude corresponds to a cluster 
mass of ~ 2 X 10^ Mq. Vast efforts have been undertaken 
from the theoretical point of view in order to explain the 
underlying universal mass function. The many publica- 
tions on this topic can be separated into those trying 
to identify some particular initial condition that selects 
a certain mass scale for star formation (e.g., Peebles & 
Dicke 1968, Fall & Rees 1985, West 1993), and those 
looking for a destruction mechanism that selects clus- 
ters in a particular mass range starting from an initially 
wide mass spectrum (e.g., Fall & Rees 1977, Gnedin & 
Ostriker 1997, Prieto & Gnedin 2008) 

At the high-mass end (i.e. nigc < /i) the mass function 
of globular clusters resembles very closely the mass func- 
tion of young clusters and molecular clouds in the Milky 
Way and other nearby galaxies (see e.g. Harris & Pudritz 
1994, Elmegreen & Efremov 1997, Gieles et al. 2006). 
On the other hand, neither young clusters nor molecu- 
lar clouds show a turnover on their mass distributions, 
but they keep rising monotonically following a power-law 
to lower masses. Fall & Zhang (2001) used simple ana- 
lytical models (including evaporation by two-body relax- 
ation, gravitational shocks and mass loss by stellar evo- 
lution) to study the evolution of the GC mass function. 
They showed that, for a wide variety of initial conditions, 
an initial power-law mass function develops a turnover 
that, after 12 Gyr, is remarkably close to the observed 
turnover of the GCLF. Vesperini (2000, 2001) reaches 
a similar conclusion, but finds that a log-normal mass 
function provides a better fit to the data. Fainter than 
the turnover, the evolution would be dominated by two- 
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body relaxation, and the mass function would end up 
having a constant number of GCs per unit mass, reflect- 
ing the fact that the masses of tidally limited clusters 
are assumed to decrease linearly with time until they 
are destroyed (other authors propose different mass-loss 
rates, see e.g., Lamers et al. 2006). Brighter than the 
turnover, the evolution is dominated by stellar evolution 
at early times and by gravitational shocks at late times. 
Recently, McLaughlin & Fall (2008) have shown that the 
GC mass function in the Milky Way depends on cluster 
half-mass density (i.e. the mean density within a radius 
containing half the total mass of the GC), in the sense 
that the turnover mass increases with half-mass density, 
while the width of the GC mass function decreases. But 
while there is currently a fairly good understanding of 
the dynamical processes that shape the GCLF, many 
details arc still missing. In particular none of the the- 
ories proposed has been entirely successful on addressing 
the question of how the turnover magnitude can remain 
constant regardless of environmental properties and the 
mass of the host galaxy. 

The use of deep HST data during the last years has re- 
sulted in high quality GCLF data, reaching 2 magnitudes 
beyond the turnover at the distance of the Virgo cluster 
(~16.5 Mpc, Mei et al. 2007). The use of these deeper ob- 
servations has recently uncovered a strong correlation be- 
tween the GCLF dispersion and the absolute magnitude 
of the parent galaxy (Jordan et al. 2006, 2007b), demon- 
strating the non-universality of this parameter and, as a 
consequence, of the GCLF as a whole. Here we present 
a study of the GCLF of 132 early type galaxies aimed 
to perform a precise test of the GCLF as a distance in- 
dicator by comparing the relative distance between the 
Virgo and Fornax clusters derived using the GCLF to 
the one derived using an analysis of surface brightness 
fluctuations (SBF, Tonry & Schneider 1988) based on 
the same data (Blakeslce et al. 2009). Previous pa- 
pers in the this series have presented an introduction to 
the survey (Jordan et al. 2007a), the properties of the 
central surface brightness proflles of early-type galax- 
ies (Cote et al. 2007) and a catalog of SBF distances 
and a precise measurement of the Virgo-Fornax distance 
(Blakeslee et al. 2009). 

The organization of this paper is as follows. In §2 wc 
present a description of the observations and data re- 
duction procedures. In §3 we describe the GCLF model 
fitting, and in §4 we compare the properties of the fits to 
previous results regarding the dispersion of the GCLF. 
Section §5 is focused on determining how universal the 
value of the turnover magnitude is, while in §6 we look for 
a better understanding of the external parameters that 
might affect this value. Finally, in §7 we summarize our 
results and the main conclusions of this paper. 

2. DATA AND GCLF INGREDIENTS 

Each one of the 132 galaxies included in this study was 

observed with the Advanced Camera for Surveys (ACS) 
during a single Hubble Space Telescope (HST) orbit, as 
part of the ACS Virgo Cluster Survey (ACSVCS) and 
the ACS Fornax Cluster Survey (ACSFCS). The goals 
and main observational features of these two surveys are 
extensively discussed in Cote et al. (2004) and Jordan 
et al. (2007a), respectively. We refer the interested reader 
to these publications for further details. 



The surveys targeted a total of 100 galaxies in the 
Virgo cluster and 43 galaxies in Fornax, and included 
observations in the F475W (« Sloan g) and F850LP (« 
Sloan z) passbands, with exposure times of ~750s and 
-^12108 respectively. In what follows we will refer to the 
F475W filter as "5" and to F850LP as due to their 
close proximity to the corresponding Sloan passbands. 

Jordan et al. (2004) describes the pipeline implemented 
to automate the reduction procedure and analysis of 
all images in both surveys. The final output from this 
pipeline is a preliminary catalog of GC candidates and 
expected contamination per galaxy, including photomet- 
ric and morphological properties, that are later used to 
evaluate the probability poc that a given object is a GC 
(see Jordan et al. 2009 for details). For the purposes 
of this study, and as defined on previous ACSVCS and 
ACSFCS papers, we constructed the GC candidate sam- 
ples by selecting all som'ces that have poc >0.5. 

Our catalog of GC candidates in a given galaxy differs 
from the intrinsic GC population due to two efl^ects: the 
existence of contamination in the sample and the level of 
completeness of the observations. 

In order to quantify the average number of contami- 
nants per field of view we have used archival ACS imaging 
of 17 blank- high latitude fields that have been observed 
in both the g and z bands, to the same or deeper depth 
than our images. These control fields were processed us- 
ing the same pipeline implemented for the science data, 
and were then used to build customized control fields, as 
if a given galaxy was in front of it (the details of this 
process are explained in Peng et al. 2006, were also a full 
list of the control fields used is available). For each of 
our target galaxies, the result is a catalog containing 17 
different estimates of the expected foreground and back- 
ground contamination. These are later used to obtain 
an average estimate of the contamination in the field of 
view of a given galaxy. 

The completeness function needs to be built consider- 
ing four parameters: the magnitude of the source (m), 
its size as measured by the projected half light radius 
(rh), its color {{g — z)o), and the surface brightness of 
the local background over which the object lies (/{,). The 
completeness function /(m, r^, {g — z)o, If,) was obtained 
by performing simulations that added model GCs of dif- 
ferent sizes (r/j = (1,3,6,10) pc), colors {{g — z)o = 
(0.7, 1.1, 1.5, 1.9) mag), and with King (1966) concentra- 
tion parameter of c = 1.5 to the images. Although the 
effect of the color of the clusters has not been considered 
in previous publications (e.g. Peng et al. 2006, Jordan 
et al. 2007b), we have now established that it also has a 
small but measurable effect over the expected complete- 
ness. Overall, roughly 6 million fake GCs were added for 
the completeness tests for each color, with equal fractions 
at each of the four sizes and avoiding physical overlaps 
with sources already present. These images were then 
reduced through exactly the same procedure used with 
the science data. The final output of the process is a 
four dimensional table that is used to evaluate / given 
an arbitrary set of (m, rh, (g — z)o, h)- The random un- 
certainty in the mean completeness curve is essentially 
zero, so the completeness limits at 90% and 50% are ro- 
bust and can be determined with negligible error for a 
given population of objects. 

This paper focuses on the study of the 89 early-type 
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Fig. 1. — Left: GCLF histogram for VCC1226 as presented in 
Jordan et al. (2007b). The hnes show the best-fit model (solid 
black curve), the intrinsic Gaussian component (dashed curve), 
the Gaussian component multiplied by the expected completeness 
(dotted curve), and a kernel density estimate of the expected con- 
tamination (solid gray curve). Right: The same as shown in the 
left hand side, but now using the corrected completeness function. 

galaxies discussed by Jordan et al. (2007b) and all 43 
galaxies of the ACSFCS. Our analysis is restricted to 
those galaxies that have more than five GC candidates 
and for which we were able to usefully constrain the 
GCLF parameters. These restrictions exclude 11 galax- 
ies in the Virgo sample but none in Fornax. 

3. GCLF MODEL FITTING 

Given the observational information previously de- 
scribed we aim to recover the parameters of the intrinsic 
luminosity function of the GCs in a galaxy. We used 
a maximum likelihood approach similar to the one de- 
scribed by Seeker & Harris (1993). According to this 
formalism, and as detailed in Jordan et al. (2007b), we 
describe the intrinsic GCLF by some function G{m \ Q), 
with Q being the set of model parameters to be fitted, 
and we assume that the uncertainties on magnitude mea- 
surements ern have a Gaussian distribution. In absence 
of contamination, the probability of observing a GC with 
a given effective radius Rh and apparent magnitude m 
against a galaxy background h would be: 



GT{m\Q,Rh,Ib, Em) = A[ h{m\em) (g) G(m \ 6) ]f{m,Rh,Ib), 

(1) 

where h{m \ £„) = (27re^) ^/^exp(-rn^/2e^), is the 
magnitude error distribution, which is convolved with 
the intrinsic GCLF G(m | 8). The normalization factor 
^ is a function of the GCLF parameters 9 and the GC 
properties Rh-, h, and e„i, and it is set by requiring that 
Gt integrates to unity over the whole magnitude range 
covered by the observations. 

In practice a fraction B of the sources classified as 
GC candidates in a galaxy are contaminants, so that 
the probability of observing a GC with parameters 
(m, /fc. Cm) is reduced by a factor (1 — B) and the dis- 
tribution that accounts for all the observed objects has 
to include the contaminants luminosity function b{m). 
Thus, the likelihood of observing a total number of N 
objects with magnitudes rrii and properties (Rh, h, ^m) 
is 



N 



£(9, B) = [][(l-Z?)GT(m, 1 9, Rha, h,^,ern,^)+Bb{m,)l 

(2) 
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Fig. 2. — Difference in turnover magnitude produced by using 
the completeness function presented by Jordan et al. 2007b and 
the one we are using here (A/i = fioiii — /J-new) in the g (top) 
and z (bottom) bands, vs. the _B-band apparent magnitude of the 
parents galaxy. 



Jordan et al. (2007b) have made a detailed description 
of several parametrization of the GCLF and their various 
advantages and drawbacks. Here we focus on the study 
of the Gaussian representation, because of its historic 
use in the study of the GCLF as a distance indicator. 
It is worth noticing that other parametrization such a 
function have also been successfully used for this purpose 
(Seeker 1992, Kissler et al. 1994). For the case of a Gaus- 
sian the set of model parameters will be 9 = {/i,(T,„}, 
where /i and a are the turnover and the dispersion in a 
distribution of the form: 



dN 
dm 



1 



2tt(t 



exp 



(m — ^)^ 



2cr2 



(3) 



The coding implementation of the outlined maximum 
likelihood procedure is in practice the same used to com- 
pute the GCLF by Jordan et al. (2007b)", except that 
we are now using completeness curves customized to the 
Fornax data too. Also, during the analysis of the ACS- 
FCS data we found a coding mistake in the interpola- 
tion of the completeness curves previously used to esti- 
mate the GCLF parameters of the Virgo galaxies. The 
background information in the completeness curves was 
sometimes misread in such a way that the completeness 
level assigned to a given background brightness was lower 
than the real value. As the changes in completeness are 
more significant for brighter backgrounds, massive galax- 
ies were more affected than dwarf galaxies. Even though 
it does not have any significant effect over the main con- 
clusions of Jordan et al. (2007b), we are reporting the 
problem here because it produces a slight change in the 
turnover magnitudes of the Virgo galaxies. The mas- 
sive galaxies are the most affected, with their turnover 
magnitudes becoming roughly '^O.l mag brighter. This 
behavior can be observed in Figure 1, where we have 
plotted side-by-side the z-band GCLF fit for VCC1226 
as presented in Figure 4 of Jordan et al. (2007b), and the 
current fit implemented using the corrected completeness 
function that now also includes a color correction. In Fig- 
ure 2 we have plotted the observed change in the turnover 
magnitude (A/i = pioid~ l^new) in both bands, against the 
B-band apparent magnitude of the parent galaxy, show- 
ing that the brightest galaxies are the most evidently af- 

In §4.2 of Jordan et al. (2007b) we show using simulations that 
our fitting procedures lead to no significant biases in the recovered 
and (T for the range of GC system sizes in our sample. 
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Fig. 3. — GCLF histograms for the Virgo and Fornax sample 
galaxies. For each one of them we present the z- and g-band GCLFs 
side by side. The VCC/FCC name and _B-band magnitude of the 
galaxy are indicated in the upper left corner of the left panel, where 
we also indicate the total number of sources A'^ in each histogram 
and the bin width h used to construct it (h is calculated as de- 
scribed in the text). In addition, we show the best-fit model (solid 
black curve), the intrinsic Gaussian component (dashed curve), 
the Gaussian component multiplied by the expected completeness 
(dotted curve), and a kernel density estimate of the expected con- 
tamination in the sample (solid gray curve) . The solid black curve 
is the sum of the solid gray and dotted curves. The galaxies are 
ordered by decreasing apparent B-band total luminosity, reading 
down from the upper left corner. The parameters of the fits are 
given in Table 1 and 2. The full version of this figure is published in 
the electronic edition of the Astrophysical Journal. Sample panels 
are shown here for guidance regarding its form and content. 

fected, unlike the dwarfs whose turnover stays virtually 
unchanged. Some spread can be observed in the case of 
the intermediate-luminosity galaxies, but in all cases the 
change in /i is always lower than 0.15 mag. 



Table 1 lists the corrected values for the Gaussian 
GCLF parameters of the ACSVCS galaxies. Updated 
values for the evolved-Schechter function fits presented 
by Jordan et al. (2007b) will be presented elsewhere. The 
Gaussian parameters shown in Table 1 are the ones con- 
sidered for this publication and they should be used for 
future reference. This table includes, for all the ACSVCS 
galaxies: the i?-band apparent magnitude from Binggeli 
et al. (1985), the estimated GCLF parameters in both 
bands, the fraction of objects that are considered to be 
contaminants, and the total number of globular cluster 
candidates (including contaminants). Table 2 presents 
the equivalent information computed for the ACSFCS 
galaxies, including the B-band absolute magnitude from 
Ferguson (1989a). Figure 3 shows the z and g-band 




0.5 1 1.5 Z 0.9 1 1.1 1.2 1.3 



<g-z> 

Fig. 4. — Left; Estimate of Gaussian dispersion in the z band, 
(Tz, vs. the same quantity in the g band, ag , for the GCLFs 
of our Fornax sample. Uncertainties are \a. The line marks the 
one-to-one correspondence between these two quantities. Right: 
Difference between estimates of Gaussian means in the g and z 
bands, fig — /iz , vs. the mean color {g — z) of the GC systems 
of our sample galaxies. Uncertainties are la. The line marks the 
one-to-one correspondence between these two quantities. 

GCLF histograms of the sample galaxies, ordered by de- 
creasing apparent B-band total luminosity. The dashed 
curve corresponds to the intrinsic Gaussian component 
given by Equation 3 and the parameters in Table 2. The 
Gaussian component multiplied by the expected com- 
pleteness is represented by the dotted curve, and a kernel 
density estimate of the expected contamination in the 
sample appears as a solid gray curve. The solid black 
curve is the sum of the solid gray and dotted curves, 
and corresponds to the net distribution for which the 
likelihood in Equation 2 is maximized. The name and 
apparent B magnitude of the galaxy are indicated in the 
upper left corner of the left panel, where we also quote 
the total number of sources in each histogram and the 
bin width h. The width of the bins, used only for dis- 
play purposes here, follows the rule h = 2(IQR)iV^^/'^ , 
where (IQR) is the interquartile range of the magnitude 
distribution and is the total number of objects in each 
GC sample (Izenman, 1991). 

As a sanity check of our fitting procedure, in the left- 
hand side of Figure 4 we compare the Gaussian disper- 
sion inferred from the GCLF fit in each band, Gg vs. 
(72, including only data from the Fornax sample. In the 
right-hand side of the same figure we have plotted the 
difference between estimates of Gaussian means in the g 
and z bands (/ig — /i^), vs. the mean color (g — z) of 
the GC systems of our sample galaxies. From the very 
tight correlation between the measurements in different 
bands we conclude that the GCLF fitting procedure is 
internally consistent and also that our error estimations 
are realistic. 

4. THE cr - MB^gai RELATION 

One of the main results discussed in Jordan 
et al. (2006, 2Q07b) is the existence of a strong corre- 
lation between the dispersion of the GCLF cr, and the 
i?-band absolute magnitude of the host galaxy Ms.gai, 
with brighter galaxies showing higher dispersion values. 
Even though some suggestive evidence on this respect 
was previously presented by other authors (e.g. Kundu 
& Whitmore, 2001) the high precision and homogene- 
ity of our ACS/HST data unveiled the a — MB,gai cor- 
relation as a general trend in GC systems, which was 
later extended to still higher galaxy luminosity by Har- 
ris et al. (2009) using 5 giant elliptical galaxies in the 
Coma cluster. Figure 5 shows this correlation for all 
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the 132 galaxies in our sample in both bands, now us- 
ing the homogeneous z-band absolute magnitudes de- 
rived from the apparent magnitudes estimated by Fer- 
rarese et al. (2006) and Cote et al. (2010, in prepara- 
tion) and the corresponding distance moduli published 
by Blakeslee et al. (2009). These values were corrected 
for reddening assuming = 1.485 E{B — V) (Ferrarese 
et al. 2006) where the value of E{B — V) was taken from 
Schlegel et al. (1998). In this figure we have used differ- 
ent symbols in order to identify the galaxies according 
to their morphological classification, but no particular 
trend related to this property seems to be obvious. The 
straight lines drawn in the panels correspond to error- 
weighted linear characterizations of these trends: 



a, = (1.07 ± 0.02) - (0.10 ± 0.01)(M,,<,ai + 22) (4) 



and 



ag = (1.10 ± 0.01) - (0.10 ± 0.01)(Af,,gai + 22). (5) 

We have excluded from these fits three galaxies for 
which no z-band magnitudes are available: VCC1535, 
VCC1030, and FCC167. Ahhough shown in Figure 5, 
FCC21 (=NGC1316) is also not included on the fits be- 
cause the observed GC system in this galaxy is highly 
influenced by interaction and proximity with its satel- 
lite galaxies, and therefore our GCLF fit is not reliable. 
Unlike Jordan et al. (2006, 2007b) we have now also ex- 
cluded from the analysis four galaxies in the Virgo clus- 
ter (VCC1297, VCCil99, VCC1192 and VCC1327) and 
two galaxies in Fornax (FCC202 and FCC143) because 
their GC systems appear to be contaminated by their 
proximity to massive ellipticals. All these galaxies are 
nonetheless retained in the Tables for completeness. 

Equations 4 and 5 confirm the trend previously ob- 
served in Virgo, and with higher statistical significance, 
by including the Fornax data. This result shows by itself 
that the GCLF parameters are not universal and depend 
at least on one parameter, i.e. the luminosity of the par- 
ent galaxy, adding an additional feature that needs to be 
accounted for by theories aiming to explain the shape of 
the GCLF. 

When the data corresponding to each cluster are fitted 
independently the linear characterizations obtained are, 
in the case of Virgo: 
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Fig. 5. — Top: GCLF dispersion cr^, inferred from Gaussian fits 
to the 2;-band data, vs. galaxy M^ g^i- The dashed line corresponds 
the linear relation between CTz and M^ g^i in eq. 4. Bottom: Same 
comparison, but for the Gaussian dispersion of the g-band GCLFs, 
cFg. The dashed line represents eq. 5. In both panels the black 
symbols correspond to Virgo galaxies, and the red ones to the 
sample in the Fornax cluster. We have morphologically separated 
the galaxies into elliptical (circles), lenticular (triangles), and dwarf 
(squares) galaxies. 



a, = (1.09 ± 0.01) - (0.08 ± 0.01)(M,,g„, + 22), (6) 
Og = (1.11 ± 0.01) - (0.09 ± 0.01)(M,,3,, + 22); (7) 

and for the Fornax cluster: 



a, = (1.07 ± 0.04) - (0.13 ± 0.02)(M^,gai + 22), (8) 
ag = (1.09 ± 0.03) - (0.10 ± 0.01)(M,,g,i + 22). (9) 

This translates into a 0.05 — 0.1 mag difference in dis- 
persion at Mz ^ 22, and also shows that the linear fits 
derived from both sets of data are equivalent within the 
uncertainties. 

As discussed in Jordan et al. (2006) it is rather straight- 
forward to link this trend in luminosity dispersion with 
a similar trend in the mass distribution of GCs. It 



is well known that giant galaxies tend on average to 
have more metal-rich GC populations when compared to 
dwarfs, showing also larger dispersions in metallicity (see 
e.g. Peng et al. 2006). This result, added to the depen- 
dence of the cluster mass-to-light ratios (T) on metallic- 
ity, opens the possibility that the observed dispersion in 
the value of a might be metallicity-driven. These vari- 
ations in T have a strong dependence on wavelength. 
In bluer filters (the g-band in our case), variations of a 
factor of 2 or more in T can be observed in the typical 
metallicity range of GCs (-2 < [Fe/H] < 0). At red- 
der wavelengths this variation becomes less dramatic, as 
shown by old stellar population models (e.g. PEGASE 
population synthesis models; Fioc & Rocca-Volmerange, 
1997). In particular the expected variation in cr as a con- 
sequence of changes of T in our z-band measurements 
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should not be higher than ^ 4%, which means that the 
spread in the value of a observed in the upper panel of 
Figure 5 reflects almost entirely a trend in the mass dis- 
tribution of globular clusters. Moreover, the very similar 
values obtained for a in the z- and (/-bands shows imme- 
diately that the trend of a with Mb cannot be generated 
by metallicity-driven changes in T. 

5. A RELATIVE VIRGO-FORNAX DISTANCE ESTIMATION 

Several methods have been used in order to obtain ac- 
curate distance estimations for both the Virgo and For- 
nax clusters, a task that is in general more easily achieved 
in the case of Fornax due to its more compact nature. 
The Virgo cluster extends for over 100 deg^ in the sky, 
showing a complex and irregular structure, with galaxies 
of different morphological type showing different spatial 
and kinematic distributions. Working on these condi- 
tions the various distance estimators have reached differ- 
ent levels of accuracy (see Ferrarese et al. 2000a, 2000b). 
We will discuss now a compilation of results from the 
literature, which are also summarized in Table 3. 

The HST Key Project to measure the Hubble con- 
stant aimed at obtaining accurate distances to galaxies 
using the period-luminosity relation for Cepheid vari- 
ables (their final results are presented in Freedman et 
al. 2001). It included the identification of Cepheids be- 
longing to 6 spiral galaxies in Virgo and 2 in Fornax, that 
were used to estimate the distance to their parent galax- 
ies, and then to the corresponding clusters. This resulted 
in distance moduli of (m — M)y = 30.92 ± 0.05 mag and 
(to — M)p = 31.39 ± 0.20 mag for Virgo and Fornax re- 
spectively, which translates into a relative distance mod- 
ulus of A(to - M)= 0.47 ± 0.20 mag. 

D'Onofrio et al. (1997) derived the relative distance 
between Virgo and Fornax by applying the D„ — cr 
(Dressier et al., 1987) and the fundamental plane (Djor- 
govski & Davis, 1987) relations to a homogeneous sample 
of early-type galaxies. The two distance indicators gave 
consistent results with a relative distance modulus of 
A(to — M) = 0.45 ± 0.15 mag. These results are in close 
agreement with the value A(to — M) = 0.52 ± 0.17 mag 
later published by Kelson et al. (2000) obtained also by 
using the fundamental plane and D„ — a relations built 
from data calibrated by the aforementioned Cepheid dis- 
tances to spiral galaxies in both Virgo and Fornax. 

The planetary nebula luminosity function (PNLF) has 
also been used for measuring distances in the local Uni- 
verse. CiarduUo et al. (1998) determined a distance 
modulus of (to - M)y = 30.79 ±0.16 mag to M87, in 
good agreement with previous measurements (e.g. Ja- 
coby et al., 1990). McMillan et al. (1993) used the 
PNLF to determine the distance to three galaxies in 
Fornax, obtaining a mean distance to the cluster of 
(m- M)f = 31.14 ±0.14 mag. If we consider M87 to be 
at the center of Virgo, the corresponding relative distance 
modulus would be A(to - M) = 0.35 ± 0.21 mag. Fer- 
rarese et al. (2000a) calibrated literature measurements 
of the PNLF using Cepheids, which led them to estimate 
a relative distance modulus between Virgo and Fornax of 
A (to — M) = 0.30 ± 0.10 mag when considering the A- 
sul)cluster as indicative of the distance to Virgo. 

Earlier relative distance modulus results derived by us- 
ing the GCLF as distance indicator present some hints 
of disagreement with the other estimations discussed 



here. Even though they were working with small and 
rather heterogeneous samples, previous studies tend to 
put this value around a very low A(to — M) ~ 0.13 mag 
(e.g. Kohle et al. 1996, Blakeslee & Tonry 1996, Ferrarese 
et al. 2000a, Richtler 2003). 

One of the most reliable distance estimators when it 
comes to population 11 samples is the surface bright- 
ness fluctuations (SBF) method due to its high internal 
precision. The ACS Virgo and Fornax clusters surveys, 
among whose aims is studying GC properties and mea- 
suring surface brightness fluctuations, provide us with 
the ideal data for comparing the properties of the GCLF 
as a distance estimator with SBF results. We will discuss 
these results separately in the next session. 

5.1. SBF distances 

The method of SBF was first introduced by Tonry 
& Schneider (1988), and uses the fluctuations produced 
in each pixel of an image by the Poissonian distribu- 
tion of unresolved stars in a galaxy in order to estimate 
the distance to the object. The amplitude of those sur- 
face brightness fluctuations normalized to the underlying 
mean galaxy luminosity are inversely proportional to dis- 
tance and can therefore be used as a distance indicator 
(see Blakeslee et al. 1999 for a review). 

The distances to the Virgo galaxies included in the 
ACSVCS have been measured using the SBF method. 
Mei et al. (2005a) describes the reduction procedure used 
for the surface brightness analysis of the ACSVCS data, 
and Mei et al. (2005b) presents the calibration for giant 
and dwarf early-type galaxies. Finally, Mei et al. (2007) 
introduces the distance catalog for a total 84 galaxies 
(50 giants and 34 dwarf) for which the SBF method was 
successfully implemented, delivering at the same time 
a three dimensional map of the structure of the Virgo 
cluster. These distance values were later updated and 
the measurements extended to include the 43 early-type 
galaxies of the ACSFCS in Blakeslee et al. (2009). In our 
analysis we will use the consistent set of Virgo and For- 
nax distances presented by the later publication. When 
no SBF distance is available for one of our sample galax- 
ies, we assume it is located at the mean Virgo distance 
((to - M) = 31.09 mag) adopted by Mei et al. (2007). 
This estimate is based on ground-based 1-band SBF mea- 
surements calibrated against Cepheids distances (Tonry 
et al. 2000, Freedman et al. 2001). 

From their SBF measurements Blakeslee et al. (2009) 
derives a relative Virgo-Fornax distance modulus 
of A(to - M) = (0.42 ± 0.03) mag, which locates 
the Fornax cluster at a distance of di? = 20 Mpc 
((to — M)f = 31.5 mag). This value is in good agree- 
ment with the relative distance moduli derived from 
the other distance estimators discussed above and 
summarized in Table 3, but it is significantly more 
precise. 

5.2. as Distance Indicator 

One of the main problems in understanding the prop- 
erties of the turnover of the GCLF as distance indicator 
is the lack of homogeneity in the data. The most com- 
prehensive compilation of recent data (mainly HST ob- 
servations) was presented by Richtler (2003), including 
a total of 102 turnover magnitudes coming from at least 
8 different publications. This inhomogeneity introduces 
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Fig. 6. — GCLF turnover magnitude, /iz, inferred from Gaussian 
fits to the 2-band data, vs. 2-band galaxy absolute magnitude 
M^^gah for all the galaxies in the ACS Virgo (black symbols) and 
Fornax (red symbols) cluster surveys. The lines represent a si- 
multaneous error-weighted linear fit performed over both samples, 
that corresponds to fi^ = (23.51 ± 0.11) + (0.04 ± 0.01)M^^gai, 
plus and offset of A(m — M) = 0.20 ± 0.04 mag for the galaxies 
in Fornax. The sample has been morphologically separated into 
elliptical (circles), lenticular (triangles), and dwarf (squares) 
galaxies. 



a major source of uncertainty in the analysis, as one has 
to rely on each author's results irrespective of the fact 
that they might not be using the same procedure to re- 
duce the data, the observations might not be on the same 
photometric band, and they might not even be using the 
same analytic form to fit the GCLF. 



The data we are presenting here are the largest and 
most homogeneous set of GCLF fits available to date. 
Our photometry is also deep enough to cover the GCLF 
at least 2 magnitudes past the turnover, therefore we 
are able to obtain reliable estimates of this parameter. 
In Figure 6 we have plotted the GCLF turnover magni- 
tude against the z-band absolute magnitude of the par- 
ent galaxy. The lines show the best linear fit p{Mz^gai) 
to each cluster's data, derived by minimizing the value 
of calculated as: 



i=l 



^^.,v-Pm,gal,v) 



A)-p(Mfw.-A) 



(10) 



where the two sums are over the ny and np galaxies 
in Virgo and Fornax respectively, S{fj,z) is the estimated 
error in /i^, and the offset A = A(m — M) corresponds 
to the relative distance modulus. In this equation each 
one of the Mz^gai components was estimated as: 



Igal = ml^gal - - M)v - A', 



(11) 



where (m — Af)v'=31.09 mag is the assumed mean dis- 
tance modulus to the Virgo cluster. The four galaxies be- 
longing to the W cloud in our Virgo sample (VCC538, 
VCC571, VCC575, VCC731 and VCC 1025) were ex- 
cluded from all our distance estimation fits as they are 
know to be located much further {D ^ 23 Mpc) than the 
mean Virgo distance. 

We have found that the best fit model for Eq. 10 cor- 
responds to = (23.51 ± 0.11) + (0.04 ± 0.01)M^,ga; 
for a value of A(m — M)= 0.20 ± 0.04 mag, where the 
error was estimated using bootstrap resampling of the 
data. This relative distance modulus represents a factor 
~2 difference with the results coming from most of the 
distance estimators previously described, and in particu- 
lar it is ~0.22 mag lower than the A(m — M)= 0.42 mag 
derived by using the SBF method with the same data. 
It is important to stress that this discrepancy cannot be 
attributed to the data itself, because we are now using a 
large sample of highly homogeneous data. Also the fact 
that the z-band absolute magnitudes of the galaxies in 
both samples were derived from equivalent observations 
and performing essentially the same analysis, minimizes 
the amount of possible biases. 

On the other hand, we are aiming to establish the level 
of precision at which /it might be useful as a distance 
indicator and therefore it seems natural to calibrate it 
against a parameter that is distance independent, which 
is not the case for Mz^gai- The GCLF dispersion, a, ap- 
pears like a good choice due to the already established 
correlation between a and Mz^gai- In Figure 7 we have 
plotted /i^ against cr^ for the complete sample in Virgo 
(black) and Fornax (red), separating the galaxies by mor- 
phological type. A minimization equivalent to Equa- 
tion 10 was also performed in this case, obtaining as the 
best fit model: ^lz = (22.99 ± 0.04) - (0.23 ± 0.04)cr^. In 
this case the offset between both samples corresponds to 
A(m — M)— 0.21 ± 0.04 mag, where the error was es- 
timated performing a bootstrap resampling of the data. 
This independent fit delivers a relative distance modulus 
that is consistent with the previously derived value. 

The observed difference between SBF and GCLF 
distances has already been reported by Richtler 
(2003), attributing this phenomena to the presence of 
intermediate- age GCs, which might contaminate the 
sample. Our sample is made up exclusively of early type 
galaxies, which are old stellar systems where the pres- 
ence of intermediate- age clusters is rarely observed (al- 
though some cases have been reported in the literature, 
see e.g. Goudfrooij et al. 2001 for the case of NGC1316 = 
FCC21, and Puzia et al. 2002 for NGC4365 = VCC731), 
so it is unlikely that this is the reason of the observed dis- 
crepancy. Ferrarese et al. (2000a) have consistently re- 
ported discrepancies between their GCLF estimated dis- 
tances and those obtained from other estimators (partic- 
ularly SBF and PNLF). They found the GCLF turnover 
in Fornax to be a full 0.5mag brighter than the value ob- 
served in Virgo. The internal errors in the GCLF mea- 
surements and the expected uncertainty due to cluster 
depth effects were not found to be enough to explain the 
scatter in their observations, suggesting the existence of 
a second parameter driving the GCLF turnover magni- 
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Fig. 7. — GCLF turnover magnitude, /iz, vs. GCLF disper- 
sion, (Tz, both inferred from Gaussian fits to the z-band data, 
for all the galaxies in the ACS Virgo (black symbols) and For- 
nax (red symbols) cluster surveys. The sample has been mor- 
phologically separated into elliptical (circles), lenticular (trian- 
gles), and dwarf (squares) galaxies. The lines represent a si- 
multaneous error-weighted linear fit to the data corresponding 
to fiz = (22.99 ± 0.04) - (0.23 ± 0.04)(Tz, with an offset of Afi^ = 
0.21 ± 0.04 for the Fornax data. 

tude. 

One obvious way to explain the observed discrepancy 
between GCLF and SBF measurements would be a mean 
age difference for the Virgo and Fornax cluster galax- 
ies (i.e., the Fornax cluster galaxies might be younger 
by some amount, leading to a brighter turnover). The 
key question, then, is determining the age difference that 
would be needed to explain the observed ~0.23 mag dif- 
ference. According to the Bruzual & Chariot (2003) 
models, for a metallicity oi Z — 0.004 and a Salpeter 
(1955) initial mass function, the observed offset would 
be consistent with an age of roughly 9 Gyr for the For- 
nax cluster when arbitrarily assuming an age of 12 Gyr 
for Virgo. This age difference would also translate into 
slightly bluer mean colors for the Fornax GCs, which 
should be on average ~0.04 mag bluer than their Virgo 
counterparts at a fixed galaxy mass. Performing a linear 
fit to the GCs mean color {g — z) vs. M^^gai correla- 
tion of our data we found that both clusters could follow 
the same trend but including an offset of 0.022 ± 0.015 
mag to redder colors in the case of Virgo. Although al- 
most consistent with zero, this value is also consistent 
with the expected color discrepancy given by the neces- 
sary age difference. The SBF technique, in which the 
fiuctuations are calibrated against a measure of the stel- 
lar populations (i.e., color), would have this difference, if 
real, accounted for. 

5.3. The Observed Dispersion on the Value of fiz 

A relatively large scatter can be observed in the 
turnover magnitude values displayed in Figure 6. In this 
subsection we want to address the question of how much 
of this dispersion is intrinsic to the sample and how much 
is the result of observational effects. The histograms in 
Figure 8 give a better illustration of this scatter, where 
we have plotted the distribution of magnitudes around 
the mean turnover magnitude of each sample, estimated 
though a 3-sigma clipping algorithm. Subtracting the 
mean turnover magnitudes of both samples we obtain: 
{p,p — jly) — 0.15 mag and [jip — p,y) — 0.17 mag, 
which delivers a first-order estimate of the relative Virgo- 




FlG. 8. — Spread in magnitude around the weighted mean 
turnover magnitude, for Virgo (top) and Fornax (bottom) in the g 
(left) and z (right) bands. 

Fornax distance modulus. We estimate the observed dis- 
persion on the right (z-band) panels of Figure 8 in 0.31 
mag and 0.28 mag, for Virgo and Fornax respectively, 
also using a 3-sigma clipping algorithm. We are more 
interested on studying the dispersion on the z-band be- 
cause is much less sensitive to metallicity variations than 
the g-hand. 

There are then three main factors driving the spread: 
cluster depth, measurement errors and the intrinsic scat- 
ter in the turnover magnitude. From their 3D map of 
the Virgo cluster, Mei et al. (2007) have determined that 
the back-to-front depth of the cluster measured from our 
sample of galaxies is 2.4±0.4 Mpc {±2a of the intrinsic 
distance distribution). At the Virgo distance this trans- 
lates into a dispersion due to line of sight effects of ~0.075 
mag. For the Fornax cluster, Blakeslee et al. (2009) esti- 
mated a depth of 2.0^0 6 Mpc (±2ct of the distance dis- 
tribution in the line-of-sight), equivalent to a dispersion 
of ^-^0.05 mag. Therefore for both clusters, the observed 
dispersion is significantly higher than the one expected 
from the cluster depth only. 

Given the observational errors and the known depths 
of the two clusters, we would like to determine whether 
there is any intrinsic dispersion in the value of fi. In order 
to do that we simulated a distribution of N galaxies (with 
N being 89 and 43 for Virgo and Fornax respectively) 
with roughly the same intrinsic turnover magnitude (we 
included a slight trend in luminosity derived from the 
lower panel of Figure 9), and we assigned them a random 
distance by using a Gaussian depth distribution with ap- 
propriate width (0.075 mag for Virgo and 0.05 mag for 
Fornax). An additional random error was added to this 
distribution based on the observed uncertainties of our 
samples. The final distribution of magnitudes was then 
used to measure the dispersion of the simulated sample 
also by using a 3-sigma clipping algorithm. This proce- 
dure was iterated 10000 times for each sample, delivering 
a mean expected dispersion in the value of fi of 0.22 mag 
for Virgo, and 0.23 mag for Fornax. These values are 
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lower than the dispersion measured in our samples, so we 
added an additional intrinsic dispersion term to the sim- 
ulations until the observed dispersion was reached. This 
difference allows for an additional dispersion of 0.21 mag 
in the case of Virgo and 0.15 mag for the Fornax cluster, 
which can not be accounted by the cluster depth or the 
observational errors alone, and therefore corresponds to 
an intrinsic dispersion in the value of /i. 

The z-band histograms shown in Figure 8 are not sym- 
metric around zero, a higher dispersion can be observed 
for positive values of {fiz — Mz)- This is consistent with 
the fact that the GCLF parameters will always be more 
precisely determined for galaxies with larger GC systems 
and they dominate the estimation of an error-weighted 
mean. As we will discuss in §6 low luminosity galaxies 
tend to show fainter turnover magnitudes and will be 
therefore located on the positive side in Figure 8, which 
combined with the larger uncertainty on the determina- 
tion of /Lto in these systems is responsible for the larger 
scatter for the positive values of {fiz — P-z)- We stress 
that, as mentioned above, in the simulations done to es- 
timate the intrinsic dispersion this slight trend of /i with 
Mz^gai is taken into account. 

6. THE UNIVERSALITY OF Mto 

The use of the GCLF as a distance indicator is based 
on the assumption of a universal value of MtO: which has 
indeed been shown to be fairly constant (within ±0.2 mag 
for massive galaxies) for a wide range of galaxy environ- 
ments. The precision and quantity of our observations 
allow to probe for potential dependencies of Mto on fac- 
tors such us the luminosity of the parent galaxy, Hubble 
type, mean color of the GC system, and environment, 
that might lurk in the observed first-order constancy of 
Mto- 

Probing for a dependence on Hubble type is important 
because the usual procedure is to use the Milky Way and 
M31 (both spiral galaxies) data in order to calibrate the 
GCLF in distant ellipticals. Our sample consists exclu- 
sively of early- type galaxies, so we cannot study the effect 
that the Hubble type might have on the value of Mto- 
However, we will discuss this later from the point of view 
of the metallicity, as the differences in the GCLF as func- 
tion of the Hubble type have been attributed to metallic- 
ity variations between the galaxies (Ashman et al., 1995). 
We will now address the influence of these factors on our 
observed non-universal GCLF. 

6.1. Luminosity 

The question of whether bright galaxies do have the 
same Mto a-s faint galaxies is particularly interesting to 
study now that the correlation between a and Mgai has 
been clearly established. Whitmore (1997) has claimed 
that dwarf ellipticals have values of M^q which are 
roughly 0.3 mag fainter than bright ellipticals, which was 
previously also mentioned by Durrell et al. (1996). In 
principle this should not represent a problem for the use 
of the GCLF as a distance indicator, as the method is 
mostly concentrated on massive galaxies which can be 
traced to larger distances. Jordan et al. (2006, 2007b) 
have also noticed that the turnover mass is slightly 
smaller in dwarf systems {Mb > — 18) compared to more 
massive galaxies (see also Miller &; Lotz (2007), showing 
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Fig. 9. — Top: GCLF absolute 2;-band turnover magnitude 
{M^q) derived from SBF distances (Blakeslee et al. 2009) vs. 
the absolute z-band magnitude of the parent galaxy (M^ g^i), 
for all the galaxies in the ACS Virgo (black) and Fornax (red) 
cluster surveys. The error-weighted linear fit corresponds to: 
Af|,Q = (-7.66 ± 0.18) -I- (0.04 ± 0.01)M^^gai- The sample has 
been morphologically separated into elliptical (circles), lenticular 
(triangles), and dwarf (squares) galaxies. Bottom: Weighted mean 
turnover magnitude (M^q) calculated in Imag wide bins (contin- 
uous lines) and over the whole magnitude range (dotted lines), for 
the Virgo (black) and Fornax (red) sample. The blue lines also 
show the corresponding weighted mean values but using the com- 
bined sample. 



that this might be partly accounted for by the effects of 
dynamical friction. 

We investigate a possible dependence of /i on M^^gai 
in Figure 9, which is equivalent to Figure 6 but with 
the observed turnover magnitudes now transformed to 
absolute turnover magnitudes using SBF distances (Mei 
et al. 2007, Blakeslee et al. 2009). The observed values 
of ^ are relatively homogeneous in the range of M^^gai 
covered by our observations, between M^.gai ~ —24 and 
Mz,gai ~ —17, however a tendency for dwarf galaxies to 
show slightly less luminous turnover magnitudes seems to 
be present. This tendency is characterized by the linear 
fit: = (-7.66 ±0.18) + (0.04 ±0.01)M^,ga!. The in- 

terpretation of this trend needs to be considered carefully 
because, due to their low luminosity, dwarf galaxies have 
smaller GC systems and the uncertainties on the deter- 
mination of uto are therefore higher. In order to lessen 
this problem, in the lower panel of Figure 9 we have plot- 
ted the weighted mean absolute magnitude in intervals 
of 1 magnitude compared to the weighted mean abso- 
lute magnitude calculated over the whole range of mag- 
nitudes. The lower luminosity bins tend to have mean 
magnitudes that are lower than the general mean both in 
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each cluster and in the combined sample. At the lower lu- 
minosity bin (in the range between —18 < Mz.gai < ^17) 
the weighted mean absolute magnitude is 0.18 mag lower 
than the general value of -8.51, and 0.3mag lower than 
the most luminous bin (-24 < M^^gai < -23). From Fig- 
ure 9 we can confirm then the trend suggested by Whit- 
more (1997) and reported by Jordan ct al.(2006, 2GG7b), 
and we conclude that the luminosity (i.e. mass) of the 
parent galaxy has an effect on determining the peak of 
the GCLF, with fainter (lower-mass) galaxies having a 
fainter GCLF turnover. 

Limiting the analysis to only the most massive galax- 
ies in the sample {Mz.gai < -21) we obtain an average 
turnover magnitude of = —8.53 mag with a dis- 
persion of 0.18 mag. These are the galaxies that could 
potentially be used as a distance indicator, and we can 
sec here that they would deliver an accurate distance 
modulus estimation within the cosmic scatter of ±0.2 
mag. There are nonetheless environmental dependen- 
cies that need to be considered before extending these 
findings to other systems because, as we can also ob- 
serve from Figure 9, the galaxies in the Fornax cluster 
show absolute turnover magnitudes that are systemat- 
ically brighter than the Virgo sample. We discus this 
point further in section §6.3. 

6.2. Color 

One of the most important requirements that a galaxy 
needs to fulfill in order to make feasible the use of its 
GCLF as a distance estimator is that its GC population 
must be old. The presence of an intermediate-age popu- 
lation will modify the GCLF by introducing clusters that 
will have brighter magnitudes than the older population. 

The GC color distribution of our sample of 89 galax- 
ies in Virgo was presented by Peng et al. (2006), where 
it was observed that on average galaxies at all luminosi- 
ties in the samples (—24 < M^^gai < ^17) appear to have 
bimodal or asymmetric GC color distributions. As dis- 
cussed in Villegas et al. (2010, in preparation) the use 
of stellar population models allow us to discard large 
age differences between red and blue GCs if we assume 
that the mass distribution of GCs does not have a de- 
pendence on [Fe/ H] inside a given galaxy. With only a 
few exceptions, the population of blue and red GC ap- 
pear to be coeval within errors for most of the galaxies, 
which lead us to concentrate on the problem of different 
metallicities between them. For giant ellipticals, this is 
also supported by previous observational studies (Puzia 
et al. 1999; Beasley et al. 2000, Jordan et al. 2002), al- 
though there are examples of massive galaxies that ap- 
pear to have formed GCs recently triggered by mergers 
(e.g. NGC 1316, Goudfrooij et al. 2001). 

With the goal of obtaining an improved calibration for 
the value of Mto, Ashman et al. (1995) studied the ef- 
fects of metallicity on the GCLF showing that changes in 
the mean metallicity of the cluster sample produce a shift 
on MtOi provided the mass distribution does not depend 
on [Fe/H]. According to Bruzual & Chariot (2003) mod- 
els, the expected change in z-band turnover magnitude, 
M|.Q, over the range of GC mean metallicity is < 0.02 
mag, which is utterly negligible considering the observa- 
tional errors. 

From a different point of view Figure 10 shows the 
correlation between turnover magnitude and mean GC 
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Fig. 10. — Turnover magnitude vs. mean GC color relationship 
for our Fornax cluster data in g (top) and 2:-bands (bottom). The 
lines show the error weighted fit to the data and have a slopes of 
0.13 ± 0.27 in the g-band and -1.01 ± 0.22 in the 2-band. 

color {g — z), in both g (top) and z-band (bottom) for 
all the galaxies in the Fornax sample. From this plot 
it can be observed that on average fig remains constant 
as a function of (g — z), but iiz tends to be brighter 
for redder GC systems. The interpretation of this plot 
presents a degeneracy between age and mass. If we as- 
sume that the Fornax galaxies, and by extension their 
GC systems, are all basically coeval, this trend can be 
explained by the fact that the z-band turnover better 
reflects mass (as it is only loosely dependent on metal- 
licity), and therefore this is an indication that galaxies 
with lower masses (as accounted by the mean metallic- 
ity of its GC system) might have less-massive turnover 
values, which translates into fainter /i^. In the g-band, 
and as a consequence also in the nearby V-band, this 
effect is canceled by the fact that the mass-to-light ra- 
tio gets lower for GCs in lower-mass, lower-metallicity 
galaxies. Therefore the historically "constant turnover 
magnitude" of the V-band GCLF might just be a conse- 
quence of the incidental cancellation of these two factors 
at this wavelength. 

6.3. Environment 

Even if we assume that the GC populations of galaxies 
of all morphological types are formed with the same ini- 
tial mass function irrespective of available gas mass and 
metallicity, there is still the environmental factor to play 
against the existence of a universal GCLF. The particular 
media in which the clusters are formed might affect their 
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evolution, shaping an environment-dependent GCLF. 

Based on data from groups and clusters of galaxies 
Blakeslee & Tonry (1996) found evidence that Mto be- 
comes fainter as the local density of galaxies increases. 
They used the velocity dispersion of groups of galaxies 
in the local universe as a density indicator in order to 
compare the values of Mto in different environments. 
Our data support the evidence presented by Blakeslee & 
Tonry (1996) in the sense they have also found a relative 
distance modulus that is too small compared to the SBF 
measurements. The trend of Mto changing as a function 
of environment (as accounted by velocity dispersion) is 
also followed by our data. However, it is important to 
mention that in spite of its lower velocity dispersion the 
Fornax cluster is denser than Virgo (Ferguson 1989b), 
and therefore the observed tendency seems to be more 
related to the total mass of the cluster than to local den- 
sity. 

Also, as discussed in §5.2 the observed discrepancy 
on the estimation of the relative Virgo-Fornax distance 
could be interpreted as a difference in mean-age between 
the stellar populations of these two clusters of galaxies, 
with an age difference of 3 Gyr being enough to explain 
this discrepancy. The combined use of high-resolution 
cosmological simulations and semi-analytic techniques 
(De Lucia et al. 2006) has shown that the faster evolu- 
tion of protocluster regions produces star formation his- 
tories that peak at higher redshift for early-type galaxies 
hosted by more massive halos. This effect would there- 
fore produce stellar populations in the Virgo clusters that 
are on average older when compared to stellar popula- 
tion belonging to a less massive galaxy cluster like For- 
nax. Mass estimates for Virgo vary substantially (e.g., 
(0.15- 1.5) X W^Mq\ Bohringer et al. 1994; Schindler 
et al. 1999; McLaughlin 1999; Tonry et al. 2000; Fouque 
et al. 2001), but it is clear that its total mass is nearly 
an order of magnitude higher than the mass of Fornax 
(~ 7 X IQ^^Mq, Drinkwater et al. 2001). The results 
presented by De Lucia et al. (2006) predict the expected 
mean-age difference for clusters of these masses to be 
~0.5 Gyr, a value that is too low to explain the observed 
difference in turnover magnitudes, but that is also de- 
pendant on the input parameters of the simulations. 

7. SUMMARY AND CONCLUSIONS 

We used AGS/HST data in order to study the GCLF of 
89 early-type galaxies in the Virgo cluster and 43 galaxies 
in the Fornax cluster, which constitute the most homo- 
geneous set of data used to date for this purpose. The 
GCLF of these galaxies was fitted by using a maximum 
likelihood approach to model the intrinsic Gaussian dis- 
tribution after accounting for contamination and com- 
pleteness effects. From the derived values of the turnover 
magnitude and the dispersion of the Gaussian fits we con- 
clude that: 

1. The analysis of 43 early- type galaxies belonging 
to the Fornax cluster shows that the dispersion of 
the GCLF decreases as the luminosity of the host 
galaxy decreases, confirming our previous results 
obtained with Virgo galaxies (Jordan et al. 2006, 
2007b). 

2. By using the GCLF turnover magnitude as a dis- 
tance indicator on our homogeneous data set we 



derive a relative distance modulus between the 
Virgo and the Fornax clusters of A(m — M)gclf — 
0.20 ± 0.04 mag, which is lower than the one de- 
rived using SBF measurements on the same data, 
A(m - M)sBF = 0.42 ± 0.03 mag. 

3. Setting the relative Virgo-Fornax distance as that 
given by SBF implies a difference in the value of 
(mto) in the two closest clusters of galaxies, sug- 
gesting that this quantity is influenced by the en- 
vironment in which a GC system is formed and 
evolves. These results support a previous study 
by Blakeslee & Tonry (1996), who fomid a corre- 
lation between GCLF turnover magnitude and ve- 
locity dispersion of the host cluster, in the sense 
that galaxy clusters with higher velocity disper- 
sions (higher masses) host galaxies with fainter 
turnovers in their GC systems. 

4. The discrepancy in the absolute magnitude of the 
GCLF turnovers in Virgo and Fornax can be ac- 
counted for if GC systems in the Fornax clusters 
were on average ~3 Gyrs younger than those in 
Virgo (thus making them brighter). Recent re- 
sults from high-resolution numerical simulations 
(e.g. Springel et al. 2005, De Lucia et al. 2006) sug- 
gest that stellar populations of Virgo-like galaxy 
clusters (high mass and high velocity dispersion) 
were formed mostly at higher redshift compared to 
less massive and lower-dispersion clusters like For- 
nax. This trend could therefore be at least partially 
responsible for the observed discrepancy in the ab- 
solute GCLF turnover magnitudes between both 
clusters. 

5. We have measured a total dispersion on the value 
of the turnover magnitude of 0.31 and 0.28 mag 
for Virgo and Fornax respectively. We show using 
simulations that these values can be only partially 
accounted by the dispersion produced by cluster 
depth and observational uncertainties. The addi- 
tional dispersion can be modeled by an intrinsic 
dispersion on the value of /io of 0.21 mag for the 
Virgo cluster and 0.15 mag for Fornax. 

6. The measured GCLF turnover is found to be 
systematically fainter for low luminosity galaxies, 
showing a ^0.3 mag decrease on dwarf systems, 
although we suffer from large uncertainties in that 
galaxy luminosity regime. The luminosity (i.e. ^ 
mass) of the parent galaxy seems to play an im- 
portant role on shaping the final form of the lumi- 
nosity distribution. This might be at least partly 
accounted for by the effects of dynamical friction if 
all other processes that contribute on shaping the 
mass function (two-body relaxation, tidal shocks, 
etc.) were to lead to a roughly constant Mto 
(Jordan et al. 2007b). 

7. Overall we find that GCLF parameters vary contin- 
uously and systematically as a function of galaxy 
luminosity (i.e. mass). The correlations we present 
here show no evidence for a dichotomy between gi- 
ant and dwarf early- type galaxies at M^ ~ —19.5 
{Mb ^ —18) in terms of their GC systems. This is 
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consistent with results presented in several reeent 
studies (e.g. Graham & Guzman, 2003; Gavazzi 
et al. 2005, Cote et al., 2006), and is at odds with 
earlier claims by Kormendy (1985). 
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TABLE 1 

Gaussian GCLF parameters for all ACSVCS galaxies. 
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22, 


,621 


± 


0, 


111 


1, 


,048 


± 


0, 


,085 


0, 


,097 


138 


vcc 


1720 


12, 


,29 


23, 


,670 


± 


0. 


127 


0, 


.797 


± 


0, 


,102 


22, 


,613 


± 


0, 


143 


0, 


,870 


± 


0, 


,115 


0, 


,141 


71 


vcc 


355 


12, 


,41 


24, 


,504 


± 


0. 


,279 


1, 


.208 


± 


0, 


.207 


23, 


,316 


± 


0. 


,206 


1, 


.027 


± 


0, 


.158 


0, 


.167 


62 


vcc 


1619 


12, 


,50 


24, 


,261 


± 


0. 


,219 


1, 


.074 


± 


0, 


.161 


23, 


,166 


± 


0. 


,234 


1, 


.082 


± 


0, 


.171 


0, 


.165 


66 


vcc 


1883 


12, 


,57 


24, 


,125 


± 


0. 


,187 


1, 


.135 


± 


0, 


.148 


22, 


,996 


± 


0. 


,166 


1, 


.064 


± 


0, 


.136 


0, 


.124 


83 


vcc 


1242 


12, 


,60 


23, 


,731 


± 


0. 


,113 





.927 


± 


0, 


.088 


22, 


,636 


± 


0. 


,120 


0, 


.983 


± 


0, 


.093 


0, 


.105 


116 


vcc 


784 


12, 


,67 


24, 


,269 


± 


0. 


,161 





.865 


± 


0, 


.123 


23, 


,102 


± 


0. 


,159 


0, 


.806 


± 


0, 


.131 


0, 


.179 


64 


vcc 


1537 


12, 


,70 


23, 


,662 


± 


0. 


240 





.977 


± 


0, 


.183 


22, 


,750 


± 


0. 


309 


1, 


,124 


± 


0, 


.232 


0, 


.256 


45 


vcc 


778 


12, 


,72 


24, 


,073 


± 


0, 


178 


1 


.052 


± 


0, 


,139 


22, 


,972 


± 


0, 


172 


1, 


,009 




0, 


,134 


0, 


,163 


74 


vcc 


1321 


12, 


,84 


24, 


,160 


± 


0, 


,225 





.926 




0, 


,168 


23, 


,153 




0, 


,222 


0, 


,919 




0, 


,166 


0, 


,198 


50 


vcc 


828 


12, 


,84 


23, 


,804 


± 


0, 


157 


1 


.045 




0, 


,142 


22, 


,787 




0, 


131 


0, 


,895 




0, 


,113 


0, 


,143 


80 


vcc 


1250 


12, 


,91 


23, 


,583 


± 


0, 


145 





.815 


± 


0, 


,111 


22, 


,609 


± 


0, 


154 


0, 


,831 


± 


0, 


,118 


0, 


,200 


54 


vcc 


1630 


12, 


,91 


24, 


,124 


± 


0, 


,326 


1 


.283 


± 


0, 


,232 


23, 


,104 




0, 


331 


1, 


,304 




0, 


,230 


0, 


,217 


57 


vcc 


1146 


12, 


,93 


23, 


,939 


± 


0, 


141 





.970 


± 


0, 


,186 


22, 


,749 


± 


0, 


127 


0, 


,890 


± 


0, 


,124 


0, 


,148 


82 


vcc 


1025 


13, 


,06 


24, 


,251 


± 


0, 


112 





.847 




0, 


,097 


23, 


,335 


dz 


0, 


136 


0, 


,938 


± 


0, 


,110 


0, 


,143 


104 


vcc 


1303 


13, 


,10 


23, 


,681 


± 


0, 


140 





.821 


± 


0, 


,106 


22, 


,793 


± 


0, 


139 


0, 


,805 


± 


0, 


,108 


0, 


,176 


61 


vcc 


1913 


13, 


,22 


23, 


,688 


± 


0, 


113 





.724 


± 


0, 


,103 


22, 


,675 


± 


0, 


117 


0, 


,738 


± 


0, 


,102 


0, 


,181 


65 


vcc 


1327 * 


13, 


,26 


23, 


,688 


± 


0, 


121 


1 


.262 


± 


0, 


,093 


22, 


,626 


± 


0, 


115 


1, 


,212 


± 


0, 


,088 


0, 


,081 


173 


vcc 


1125 


13, 


,30 


23, 


,667 


± 


0, 


127 


0, 


.781 


± 


0, 


,109 


22, 


,645 


± 


0. 


136 


0, 


,791 


it 


0, 


,109 


0, 


,179 


62 


vcc 


1475 


13, 


,36 


24, 


,073 


± 


0. 


,141 





.990 


± 


0, 


.107 


23, 


,199 


± 


0. 


,178 


1, 


.101 


± 


0, 


.133 


0, 


.137 


86 


vcc 


1178 


13, 


,37 


23, 


,609 


± 


0. 


,134 





.997 


± 


0, 


.102 


22, 


,562 


± 


0. 


,123 


0, 


.949 


± 


0, 


.090 


0, 


.124 


90 


vcc 


1283 


13, 


,45 


24, 


,049 


± 


0. 


,152 





.894 


± 


0, 


.120 


23, 


,023 


± 


0. 


,167 


0, 


.932 


± 


0, 


.129 


0, 


.170 


66 


vcc 


1261 


13, 


,56 


23, 


,962 


± 


0. 


,275 


1, 


.133 


± 


0, 


.208 


23, 


,004 


± 


0. 


,327 


1, 


.243 


± 


0, 


.238 


0, 


.217 


46 


vcc 


698 


13, 


,60 


23, 


,793 


± 


0. 


,090 





.843 




0, 


.066 


22, 


,777 


± 


0. 


,085 


0, 


.810 


± 


0, 


.062 


0, 


.105 


119 


vcc 


1422 


13, 


,64 


23, 


,625 


± 


0. 


,169 





.656 


± 


0, 


.130 


22, 


,521 


± 


0. 


,168 


0, 


.651 


± 


0, 


.127 


0, 


.258 


37 


vcc 


2048 


13, 


,81 


23, 


,450 


± 


0, 


,324 





.969 


± 


0, 


,217 


22, 


,420 


± 


0, 


,282 


0, 


,881 


± 


0, 


,194 


0, 


,420 


22 


vcc 


1871 


13, 


,86 


23, 


,520 


± 


0, 


,608 


1 


.181 


± 


0, 


,455 


22, 


,512 


± 


0, 


,604 


1, 


,154 


± 


0, 


,480 


0, 


,516 


18 


vcc 


9 


13, 


,93 


23, 


,940 


± 


0, 


,391 


1 


.086 


± 


0, 


,305 


22, 


,830 


± 


0, 


,260 


0, 


,894 


± 


0, 


,196 


0, 


,246 


34 


vcc 


575 


14, 


,14 


24, 


,847 


± 


0, 


,271 





.665 


± 


0, 


,281 


23, 


,833 


± 


0, 


130 


0, 


,333 


± 


0, 


,184 


0, 


,386 


27 


vcc 


1910 


14, 


,17 


23, 


,758 


± 


0, 


,208 


1 


.175 


± 


0, 


,161 


22, 


,630 


± 


0, 


,209 


1, 


,135 


± 


0, 


,176 


0, 


,180 


60 


vcc 


1049 


14, 


,20 


24, 


,052 


± 


0, 


,257 





.550 




0, 


,197 


23, 


,106 


± 


0, 


,396 


0, 


,634 




0, 


,268 


0, 


,487 


18 


vcc 


856 


14, 


,25 


23, 


,792 


± 


0, 


185 





.887 


± 


0, 


,156 


22, 


,768 


± 


0, 


164 


0, 


,862 


± 


0, 


,127 


0, 


,211 


50 


vcc 


140 


14, 


,30 


23, 


,992 


± 


0, 


,245 





.790 


± 


0, 


,197 


22, 


,979 


± 


0, 


,249 


0, 


,822 


± 


0, 


,182 


0, 


,329 


29 


vcc 


1355 


14, 


,31 


24, 


,554 


± 


0, 


,776 


1 


.273 


± 


0, 


,541 


23, 


,682 


± 


0, 


,732 


1, 


,161 


± 


0, 


,530 


0, 


,471 


20 


vcc 


1087 


14, 


,31 


23, 


,732 


± 


0, 


134 





.926 


± 


0, 


,101 


22, 


,713 


± 


0, 


133 


0, 


,898 


± 


0, 


,112 


0, 


,162 


68 


vcc 


1297 * 


14, 


,33 


23, 


,403 


d= 


0, 


109 


1 


.141 


± 


0, 


,082 


22, 


,299 


± 


0, 


105 


1, 


,084 




0, 


,080 


0, 


,092 


152 


vcc 


1861 


14, 


,37 


23, 


,608 


± 


0, 


,222 


1 


.015 




0, 


,185 


22, 


,572 




0, 


,206 


0, 


,937 


± 


0, 


,164 


0, 


,234 


49 


vcc 


543 


14, 


,39 


23, 


,854 


± 


0, 


196 


0, 


.692 




0, 


,139 


22, 


,792 


± 


0, 


184 


0, 


,635 


it 


0, 


,127 


0, 


,330 


28 


vcc 


1431 


14, 


,51 


24, 


,092 


± 


0. 


,171 


1, 


.054 


± 


0, 


.128 


23, 


,054 


± 


0. 


,188 


1, 


.082 


it 


0, 


.140 


0, 


.158 


71 


vcc 


1528 


14, 


,51 


23, 


,550 


± 


0. 


,137 





.720 


± 


0, 


.105 


22, 


,609 


± 


0. 


,129 


0, 


.697 


it 


0, 


.097 


0, 


.222 


49 


vcc 


1695 


14, 


,53 


24, 


,416 


± 


0. 


,401 





.962 


± 


0, 


.289 


23, 


,480 


± 


0. 


,517 


1, 


.103 


it 


0, 


.357 


0, 


.380 


22 


vcc 


1833 


14, 


,54 


24, 


,091 


± 


0. 


,223 





.695 


± 


0, 


.159 


22, 


,954 


± 


0. 


,147 


0, 


.500 


it 


0, 


.110 


0, 


.332 


28 


vcc 


437 


14, 


,54 


23, 


,933 


± 


0. 


,162 





.783 


± 


0, 


.134 


23, 


,056 


± 


0. 


,167 


0, 


.845 


it 


0, 


.131 


0, 


.229 


50 


vcc 


2019 


14, 


,55 


23, 


,551 


± 


0. 


220 





.873 


± 


0, 


,200 


22, 


,619 


± 


0. 


,225 


0, 


.860 


it 


0, 


.193 


0, 


.303 


34 


vcc 


200 


14, 


,69 


24, 


,459 


d= 


0, 


,221 





.680 


± 


0, 


,144 


23, 


,582 


± 


0, 


331 


0, 


,834 


it 


0, 


,221 


0, 


,381 


25 


vcc 


571 


14, 


,74 


24, 


,392 


± 


0, 


,543 





.951 




0, 


,346 


24, 


,249 


± 


1, 


,542 


1, 


,421 


it 


0, 


,810 


0, 


,478 


17 


vcc 


21 


14, 


,75 


24, 


,073 


± 


0, 


,636 


1 


.438 


± 


0, 


,418 


22, 


,963 


± 


0, 


,559 


1, 


,276 


it 


0, 


,387 


0, 


,351 


26 


vcc 


1488 


14, 


,76 


24, 


,146 


± 


0, 


,303 





.580 




0, 


,208 


23, 


,088 


± 


0, 


,390 


0, 


,553 


it 


0, 


,262 


0, 


,471 


19 


vcc 


1499 


14, 


,94 


24, 


,562 


± 


0. 


601 


1, 


.418 


± 


0, 


,377 


23, 


,489 


± 


0. 


,530 


1, 


,295 


it 


0, 


,341 


0, 


,272 


35 



GCLF as distance indicator 
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TABLE 1 

Gaussian GCLF parameters for all ACSVCS galaxies. 



vcc 


1545 


14, 


.96 


24, 


,099 


± 


0, 


164 


0, 


,910 


± 


0, 


128 


23, 


,159 


± 


0, 


183 


0, 


,9,30 


± 


0, 


145 


0, 


,189 


63 


vcc 


1192 * 


15 


.04 


23, 


,781 


± 


0. 


,086 


1, 


,073 




0. 


,066 


22, 


,663 


± 


0. 


,085 


1, 


,052 


it 


0. 


,064 


0, 


,064 


213 


vcc 


1075 


15 


.08 


23. 


.514 


± 


0. 


,169 


0, 


,554 


± 


0. 


,119 


22, 


,522 


± 


0. 


,155 


0, 


.515 


± 


0. 


,115 


0, 


.378 


26 


vcc 


1440 


15 


.20 


24. 


.267 


± 


0. 


,237 


0, 


,895 


± 


0. 


,176 


23, 


,270 


± 


0. 


,221 


0, 


.824 


± 


0. 


,162 


0, 


.259 


38 


vcc 


230 


15 


.20 


23. 


.941 


± 


0. 


,134 


0, 


,541 


± 


0. 


,106 


23, 


,078 


± 


0. 


,139 


0, 


.578 


± 


0. 


,105 


0, 


.274 


38 


vcc 


2050 


15 


.20 


23. 


.900 


± 


0. 


,118 


0, 


,281 


± 


0. 


,089 


22, 


,963 


± 


0. 


,135 


0, 


.304 


± 


0. 


,106 


0, 


.459 


20 


vcc 


751 


15 


.30 


23, 


,525 


± 


0. 


,191 


0, 


,504 


± 


0. 


,130 


22, 


,699 


± 


0. 


206 


0, 


,509 


± 


0. 


,130 


0, 


.495 


17 


vcc 


1828 


15 


.33 


23, 


,807 


± 


0. 


210 


0, 


,702 


± 


0. 


,183 


22, 


,757 


± 


0. 


,198 


0, 


,664 


± 


0. 


,148 


0, 


.355 


27 


vcc 


1407 


15 


.49 


24, 


,397 


± 


0, 


123 


0, 


,665 


± 


0, 


,094 


23, 


,420 


± 


0, 


144 


0, 


,745 


± 


0, 


111 


0, 


,186 


60 


vcc 


1886 


15 


.49 


23, 


,034 


± 


0, 


,715 


0, 


,971 




0, 


,463 


21, 


,565 


± 


0, 


,304 


0, 


,463 


± 


0, 


,215 


0, 


,622 


14 


vcc 


1199 * 


15, 


.50 


23, 


,833 


± 


0, 


,094 


1, 


,166 




0, 


,074 


22, 


,682 


± 


0, 


,089 


1, 


,125 


± 


0, 


,070 


0, 


,060 


228 


vcc 


1539 


15, 


.68 


23, 


,813 


± 


0, 


182 


0, 


,831 


± 


0, 


168 


22, 


,820 


± 


0, 


199 


0, 


,901 




0, 


147 


0, 


,275 


43 


vcc 


1185 


15 


.68 


23, 


,843 


± 


0, 


172 


0, 


,693 


± 


0, 


116 


22, 


,910 


d= 


0, 


155 


0, 


,639 


± 


0, 


105 


0, 


,292 


33 


vcc 


1489 


15 


.89 


23, 


,977 


± 


0, 


150 


0, 


,378 


± 


0, 


129 


23, 


,157 


d= 


0, 


,279 


0, 


,484 


± 


0, 


,469 


0, 


,417 


22 


vcc 


1661 


15, 


.97 


24, 


,040 


± 


0. 


,281 


0, 


,630 


± 


0. 


,273 


23, 


,058 


± 


0. 


,285 


0, 


,614 


± 


0. 


,215 


0, 


,477 


19 



Notes — (1) Galaxy VCC number. (2) Galaxy B-band magnitude. (3) and (4) Maximum likelihood estimates of the Gaussian mean 
fi and dispersion a of the g-band GCLF. (5) and (6) Same as cols. (3) and (4), but for the z band. (7) Fraction of the sample that 
is expected to be contamination. (8) Total number N of all objects (including contaminants and uncorrected for incompleteness) with 
PGC > 0.5. (*) These galaxies were excluded from the analysis because of their close proximity to massive elliptical galaxies. 
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Villegas et al. 



TABLE 2 

Gaussian GCLF parameters for all ACSFCS galaxies. 



ID 








IJ.Z 




/9 


N 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



FCC 21 


9, 


.4 


26, 


,350 


± 


1, 


.234 


2, 


.178 




0, 


.059 


25. 


,150 


± 


0, 


.668 


2, 


.189 


± 


0, 


.060 


0, 


.011 


647 


FCC 213 


10, 


.6 


24, 


,090 


± 


0, 


.048 


1, 


.231 


± 


0, 


.038 


22. 


,802 


± 


0, 


.044 


1, 


.198 




0, 


.035 


0, 


.015 


1074 


FCC 219 


10, 


.9 


24, 


,140 


± 


0, 


.072 


1, 


.110 


± 


0, 


.058 


22. 


,940 




0, 


.072 


1, 


.112 


± 


0, 


.058 


0, 


.039 


380 


NGC 1340 


11, 


,2 


24, 


,384 


± 


0, 


.098 


1, 


.124 


± 


0, 


.074 


23. 


,468 


± 


0, 


.111 


1, 


.180 


± 


0, 


.082 


0, 


.039 


280 


FCC 167 


11, 


.3 


24, 


,023 


± 


0, 


.059 


1, 


.022 


± 


0, 


.046 


22. 


,808 


± 


0, 


.060 


1, 


.044 




0, 


.047 


0, 


.026 


424 


FCC 276 


11, 


.8 


24, 


,032 


± 


0, 


.070 


1, 


.102 




0, 


.063 


22. 


,960 


± 


0, 


.076 


1, 


.166 


± 


0, 


.061 


0, 


.040 


361 


FCC 147 


11, 


.9 


24, 


,077 


± 


0, 


.085 


1, 


.197 


± 


0, 


.067 


22. 


,894 


± 


0, 


.081 


1, 


.156 


± 


0, 


.064 


0, 


.047 


320 


IC 2006 


12, 


.2 


24, 


,076 


± 


0, 


.092 


0, 


.886 


it 


0, 


.070 


22. 


,935 


± 


0, 


.089 


0, 


.868 


± 


0, 


.067 


0, 


.085 


132 


FCC 83 


12, 


.3 


24, 


,026 


± 


0, 


.076 


1, 


.040 


± 


0, 


.058 


22. 


,906 


± 


0, 


.070 


0, 


.988 


± 


0, 


.054 


0, 


.044 


274 


FCC 184 


12, 


.3 


23, 


.956 


± 


0, 


.067 


1, 


.029 


± 


0, 


.054 


22. 


,664 


± 


0, 


.067 


1, 


.030 


± 


0, 


.054 


0, 


.042 


306 


FCC 63 


12, 


.7 


24. 


.023 


± 


0. 


.106 


1, 


.236 




0. 


.084 


22. 


.951 


± 


0, 


.108 


1, 


.233 


± 


0. 


.086 


0, 


.058 


231 


FCC 193 


12, 


.8 


23. 


.934 


± 


0. 


.161 


0. 


.822 


± 


0. 


.123 


22. 


.830 


± 


0, 


.172 


0, 


.899 




0. 


.130 


0. 


.176 


48 


FCC 153 


13, 


.0 


24. 


.066 


± 


0. 


.175 


0. 


.947 


± 


0. 


.135 


23. 


.086 


± 


0. 


.155 


0. 


.857 




0. 


.119 


0. 


.161 


60 


FCC 170 


13, 


.0 


24. 


.016 


± 


0. 


.196 


1. 


.182 




0. 


.167 


23. 


.073 


± 


0. 


.202 


1. 


.195 




0. 


.171 


0. 


.137 


71 


FCC 177 


13, 


.2 


23. 


.897 


± 


0. 


.139 


0. 


.928 




0. 


.108 


22. 


.923 




0. 


.125 


0. 


.859 




0. 


.095 


0. 


.129 


70 


FCC 47 


13, 


.3 


23. 


,993 


± 


0. 


.068 


0. 


.988 


± 


0. 


.053 


22. 


.948 


± 


0. 


.068 


0. 


.984 


± 


0. 


.054 


0. 


.044 


276 


FCC 43 


13, 


.5 


24. 


,342 


± 


0, 


.304 


1, 


.088 


± 


0, 


.238 


23. 


,261 


± 


0, 


.289 


1, 


.099 


± 


0, 


.220 


0, 


.208 


37 


FCC 190 


13, 


.5 


23. 


,934 


± 


0, 


.090 


0, 


.932 


± 


0, 


.072 


22. 


,937 


± 


0, 


.091 


0, 


.940 


± 


0, 


.073 


0, 


.071 


156 


FCC 310 


13, 


.5 


24. 


,144 


± 


0, 


.167 


0, 


.743 




0, 


.122 


23. 


,184 


± 


0, 


.169 


0, 


.736 


± 


0, 


.123 


0, 


.229 


39 


FCC 148 


13, 


.6 


23. 


,851 


± 


0, 


.134 


1, 


.012 


± 


0, 


.107 


22. 


,837 


± 


0, 


.147 


1, 


.079 


± 


0, 


.117 


0, 


.111 


86 


FCC 249 


13, 


.6 


23. 


,913 


± 


0, 


.089 


0, 


.929 


± 


0, 


.068 


22. 


,935 


± 


0, 


.091 


0, 


.939 


lb 


0, 


.070 


0, 


.078 


155 


FCC 255 


13, 


.7 


23. 


,737 


± 


0, 


.111 


0, 


.780 


± 


0, 


.089 


22. 


,714 


± 


0, 


.110 


0, 


.770 


± 


0, 


.087 


0, 


.125 


80 


FCC 277 


13, 


.8 


24. 


,244 


± 


0, 


.158 


0, 


.677 


± 


0, 


.136 


23. 


,278 


± 


0, 


.156 


0, 


.683 


± 


0, 


.121 


0, 


.199 


42 


FCC 55 


13, 


.9 


24. 


,446 


± 


0, 


.148 


0, 


.655 


± 


0, 


.111 


23. 


,441 


± 


0, 


.181 


0, 


.734 


± 


0, 


.135 


0, 


.223 


37 


FCC 152 


14, 


.1 


23. 


,485 


± 


0, 


.344 


0, 


.844 


± 


0, 


.234 


22. 


,492 


± 


0, 


.248 


0, 


.585 


± 


0, 


.179 


0, 


.456 


16 


FCC 301 


14, 


.2 


24. 


,383 


± 


0, 


.415 


1, 


.008 


± 


0, 


,289 


23. 


,605 


± 


0, 


.628 


1, 


.216 


± 


0, 


.406 


0, 


.353 


21 


FCC 335 


14, 


.2 


23. 


,026 


± 


0, 


.766 


1, 


.593 




0, 


.585 


21. 


,954 




0, 


.732 


1, 


.517 


it 


0, 


.570 


0, 


.525 


14 


FCC 143 * 


14, 


.3 


23. 


.873 


± 


0, 


.148 


0. 


.908 


± 


0. 


.114 


22. 


.929 


± 


0. 


.141 


0. 


.855 


± 


0. 


.119 


0. 


.158 


62 


FCC 95 


14, 


.6 


24. 


.154 


± 


0. 


.098 


0. 


.263 




0. 


.074 


23. 


.069 


± 


0. 


.057 


0. 


.155 


lb 


0. 


.046 


0. 


.373 


21 


FCC 136 


14, 


.8 


23. 


.968 


± 


0. 


.163 


0. 


.436 


± 


0. 


.245 


23. 


.011 


± 


0. 


.104 


0. 


.355 


lb 


0. 


.088 


0. 


.294 


25 


FCC 182 


14, 


.9 


24. 


.169 


± 


0. 


.142 


0. 


.891 


± 


0. 


.111 


23. 


.220 


± 


0. 


.178 


1. 


.008 


lb 


0. 


.136 


0. 


.145 


59 


FCC 204 


14, 


.9 


24. 


.192 


± 


0. 


.518 


0. 


.944 


± 


0. 


.392 


23. 


.599 


± 


0. 


.757 


1. 


.118 


lb 


0. 


.498 


0. 


.443 


17 


FCC 119 


15, 


.0 


25. 


.464 


± 


0. 


.946 


0. 


.972 


± 


0. 


.551 


25. 


,150 


± 


7. 


.387 


1. 


.222 


lb 


0. 


.312 


0. 


.411 


17 


FCC 26 


15, 


,0 


23. 


,208 




0, 


,143 


0, 


,441 


± 


0, 


.114 


22. 


,394 




0, 


.226 


0, 


.657 


± 


0, 


.159 


0, 


.337 


22 


FCC 90 


15, 


.0 


23. 


,953 


± 


0, 


,352 


0, 


,673 


± 


0, 


.299 


23. 


,026 


± 


0, 


.213 


0, 


.567 


± 


0, 


.211 


0, 


.370 


21 


FCC 106 


15, 


.1 


23. 


,966 


± 


1, 


.321 


1, 


.990 




0, 


.917 


23. 


,235 


± 


1, 


.879 


2, 


.223 


± 


1, 


.326 


0, 


.486 


15 


FCC 19 


15, 


.2 


24. 


,552 


± 


0, 


.459 


0, 


.750 


± 


0, 


.291 


23. 


,483 


± 


0, 


.496 


0, 


.720 


± 


0, 


.320 


0, 


.463 


16 


FCC 288 


15, 


.4 


24. 


,913 


± 


0, 


.510 


0, 


.871 


± 


0, 


.353 


24. 


,355 


± 


0, 


.706 


0, 


.893 


lb 


0, 


.682 


0, 


.426 


17 


FCC 202 * 


15, 


.3 


23. 


,996 


± 


0, 


.084 


1, 


.101 


± 


0, 


.068 


22. 


,834 


± 


0, 


.083 


1, 


.087 


lb 


0, 


.069 


0, 


.050 


232 


FCC 324 


15, 


.3 


23. 


,698 


± 


0, 


.271 


0, 


.665 


± 


0, 


.207 


22. 


,947 


± 


0, 


.274 


0, 


.758 


lb 


0, 


.192 


0, 


.384 


21 


FCC 100 


15, 


.5 


24. 


,119 


± 


0, 


.114 


0, 


.433 


± 


0, 


.098 


23. 


,366 


± 


0, 


.181 


0, 


.590 


lb 


0, 


.137 


0, 


.272 


34 


FCC 203 


15, 


.5 


24. 


,155 


± 


0, 


.361 


1, 


.184 


± 


0, 


.270 


23. 


,167 


± 


0, 


.329 


1, 


.114 


lb 


0, 


.241 


0, 


.271 


30 


FCC 303 


15, 


.5 


23. 


.479 


± 


0, 


.200 


0, 


.623 


± 


0. 


.141 


22. 


.531 


± 


0, 


.218 


0, 


.678 


lb 


0, 


.152 


0, 


.350 


22 



Notes - (1) Galaxy FCC number. (2) Galaxy _B-band magnitude. (3) and (4) Maximum likelihood estimates of the Gaussian mean 
^ and dispersion a of the (/-band GCLF. (5) and (6) Same as cols. (3) and (4), but for tiie z-band. (7) Fraction of the sample that 
is expected to be contamination. (8) Total number N of all objects (including contaminants and uncorrected for incompleteness) with 
Pgc ^ 0-5. (*) These galaxies were excluded from the analysis because of their close proximity to massive elliptical galaxies. 

TABLE 3 

Literature compilation of relative distance modulus between Virgo 
AND Fornax clusters 



Method A(m — M) Reference 



Cephcids 
Fund. Plane 

PNLF 

GCLF 



SBF 



0.47 ± 0.20 
0.45 ± 0.15 
0.52 ± 0.17 
0.35 ± 0.21 
0.30 ±0.10 
0.08 ± 0.09 
0.13 ± 0.11 
0.09 ± 0.27 
0.17 ±0.28 
0.42 ± 0.03 



1 
2 
3 
4, 5 
6 
7 
8 
6 
9 
10 



Notes - The cited references arc; (1) Frecdman et al. 2001, (2) D'Onofrio et al. 1997, (3) Kelson ct al. 2000, (4) Ciardullo et al. 1998, 
(5) McMillan et al. 1993, (6) Ferrarese et al. 2000a, (7) Kohle et al. 1996, (8) Blakeslee & Tonry 1996, (9) Richtler 2003, (10) Blakeslee 
et al. 2009. 



